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Abstract 

We investigated the production of neutralino pairs via photon-photon 
collisions in the minimal supersymmetric model(MSSM) at future lin- 
ear colliders. The numerical analysis of their production rates is carried 
out in the mSUGRA scenario. The results show that this cross section 
can reach about 18 femto barn for X1X2 P a i r production and 9 femto 
barn for X2X2 P arr production with our chosen input parameters. 
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1. Introduction 



In the supersymmetric theory [|I]], proper electroweak symmetry breaking induces the right 
properties of the lightest supersymmetric particle (LSP) to be a natural candidate for weak- 
interacting cold dark matter, which can explain many astrophysical observations]^]. The 
minimal supersymmetric standard model(MSSM)[jl| predicts that there exists an absolutely 
stable LSP. Most often the LSP in the MSSM theory is the lightest Majorana fermionic 
neutralino Xi- Therefore the production of the lightest neutralino Xi an d the second lightest 
neutralino x\ ma y be studied at present and future experiments and the detailed study of the 
neutralino sector will help us to determine which kind of the supersymmetric models really 
exists in nature. In the MSSM, the physical neutralino mass eigenstates Xi (2 = 1 ~ 4) are 
the combinations of the neutral gauginos(.B, W 3 ) and the neutral higgsinos(if°, H 2 ). m the 
two-component fermion fields ip® = — zA 3 , ipH i V'.ff )!!' w here A' is the bino and A 3 is 

the neutral wino, the neutralino mass matrix in the Lagrangian is given by 

C M = -\{^) T Y^ + h.c, (1.1) 

where the matrix Y reads 

/ Mi — m z sin 9w cos (3 mz sin 9\y sin (3 \ 

y_ M 2 m z cos 9\y cos (3 — mz cos 6w sin (3 (12) 

—m z sin 6 W cos (3 m z cos 9w cos (3 — fi ^ ' ' 

\ mz sin 6w sin (3 —m z cos 9w sin (3 —\i / 

In the above equation, the neutralino mass matrix is related to four unknown parameters, 
namely fi, M 2 , M 1 and tan (3 = v 2 /vi, ratio of the vacuum expectation values of the two Higgs 
fields, /i is the supersymmetric Higgs-boson-mass parameter and M 2 and Mi are the gaug- 
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ino mass parameters associated with the SU(2) and U(l) subgroups, respectively. In CP- 
noninvariant theories, M 1; M 2 and higgsino mass parameter \x can be complex. However, M 2 
can be real and positive without loss of generality by reparametrization of the fields. In this 
work we shall investigate neutralino pair production in framework of the MSSM while ignor- 
ing CP-violation and taking M\ and fi as being real. The matrix Y is symmetric and can be 
diagonalized by one unitary matrix iV such that Nd = N*YN + = diag(m^o, m^o, m^o,m^o) 
in order of m^o < m^o < m^o < m^o. Then the two-component mass eigenstates can be 

X? = %V°, z,j = l,...,4. (1.3) 
The proper four-component mass eigenstates are defined in terms of two-component fields 



as 



and the mass term becomes 



(< = !.•• -.4), (1.4) 



c m = - l -Y,M°ixl (1.5) 

i 

where Mj's are the diagonal elements of -/Vp. 

The future e + e~ Linear Colliders(LC) are designed to give facilities for e + e~, 77 and 
other collisions at the energy of 500 ~ 2000 GeV with a luminosity of the order 10 33 cm _2 s~ 1 
Q. For example, the proposed TESLA collider is known as a powerful tool for exploration 
the multi-handred GeV scale ||. Different options of this machine, namely, e + e~, 77, je 
and e~e~ are complementary to each other and will add essential new information to that 
obtained from the LHC. Its annual (10 7 s) 77 luminosity will be about 10 — 30 (in the 



high energy peak) with possible upgrade of luminosity by one order of magnitude. Searching 
for supersymmetric particles and determining their properties are one of the main tasks at 
future LC. In detecting the existence of neutralinos, both e + e~ and 77 collisions have clearer 
background than hardron collisions, but 77 collision at LC would have distinct advantage 
over the situation of LC operating in the e + e~ colision mode, where the resonant effects 
of Higgs bosons can be observed only at some specific center of mass energy ranges of the 
machine and the possible s-channel suppression would generally reduce the cross section 
of neutralino pair production. Because of the continuous c.m.s energy distribution of the 
colliding photons back-scattered by e~ and e + beams, the intermediate resonant effects of 
Higgs bosons could enhance the neutralino pair production rate over a rather wide colliding 
energy range at electron-positron colliders. 

So far there is no experimental evidence for neutralinos at CERN LEP2. They only set 
lower bound on the lightest neutralino mass m^o . Recent experimental reports presented that 
the mass of the lightest neutralino may be larger than 32.5 Gel^H and the lower limit of the 
chargino mass is 76.8 GeV |7J. The direct neutralino pair productions at the CERN Large 
Hadron Collider (LHC) are studied in references ||||[rC]. The production of neutralino pair 
can be produced also at the LC machine operating in both e~e + and 77 collision modes. 
Recently, it has been found that the X1X2 production rate in e + e~ collision mode can reach 



several handred femto barn 11 



In this paper we investigate the potential of direct neutralino pair production at the LC 
operating in photon-photon collision mode in framework of the MSSM with complete one- 
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loop Feynman diagrams. The numerical calculation will be illustrated in the CP-conserving 
mSUGRA scenario with five input parameters, namely mi/2, mo, Aq, fi and tan/3, where 
mi/2, mo and A are the universal gaugino mass, scalar mass at GUT scale and the trilinear 
soft breaking parameter in the superpotential respectively. From these five parameters, all 
the masses and couplings of the model are determined by evolution from the GUT scale down 



to the low electroweak scale |fT2| . The paper is organized as follows: In section 2, we give 
calculations of the neutralino pair production at the LC operating in photon-photon collision 
mode. In section 3, we discuss the numerical results of the cross sections. A short summary is 
presented in section 4. Finally the explicit expressions of form factors for s-channel diagrams 
are listed. 

2. The Calculation of e + e~ — > 77 — > XiX*j 

Neutralinos can be produced at e + e~ colliders, either in diagonal or in mixed pairs. In 
this section we calculate the processes 

e + e - _> 77 _> (i = l,2, j = 2). 

The generic Feynman diagrams contributing to the subprocess 77 — > XiX*j m t ne MSSM at 
one-loop level are depicted in Fig.l, where the exchange of incoming photons in Fig.l(a.l ~ 
6), Fig.l(b.5) and Fig.l(c.l ~ 4) are not shown. Fig.l(a.l ~ 6) are box diagrams. Fig. 1(6.1 ~ 
5) are quartic interaction diagrams. Fig.l(c.l ~ 3) represent triangle diagrams. The 
Fig. 1(6. 3 ~ 6.5), (c.l ~ 3) are also called s-channel diagrams. The Feynman diagrams 
in Fig.l include the loops of quarks(U,D), squarks (U,D), leptons(E), sleptons(-E) of three 



generations, charginos, W + boson, charged ghost particles, charged Higgs boson and gold- 
stone. Due to the Yukawa couplings strength, the contributions from the loop diagrams of 
the third generation quarks and squarks are more important than those from other diagrams. 
Here we should mention two points: (1) There is no diagram with a triangle squark(slepton) 
loop coupling with an A or G° boson in Fig.l, because the vertices of A°(G°) — q — q 
vanish JTB|]. (2) Our calculation shows that the Feynman diagrams involving quartic vertices 
7 - hP - G+(H + ) - G-(H-) and 7 - H° - G+(H + ) - G-(H~), which have similar structure 
to Fig.l(b.5), have no contribution to the cross section, therefore omitted in figure 1. The 
Z° boson intermediated s-channel diagrams similar to the diagrams of Fig.l(b.3) ~ (b.5), 
Fig.l(c.l) ~ (c.3) are not plotted in Fig.l either, since they cannot contribute to the cross 
section. For this result there are two reasons: (1) The CP-odd scalar component of the 
Z° boson does not couple to the invariant CP-even XiXj state. (2) The vector component 
of the Z° boson wave function does not couple to the initial 77 state as the result of the 
Laudau-Yang theorem. 

Since there is no tree level diagram for the neutralino pair production via photon-photon 
collisions, the calculation for this process can be simply carried out by summing all unrenor- 
malized reducible and irreducible one-loop diagrams and the results will be finite and gauge 
invariant. In this work, we perform the calculation in the 't Hooft-Feynman gauge, and take 
the CKM matrix as identity. 

We denote the reaction of neutralino pair production via photon-photon collisions as: 

l(pi^Hp2^)^Xi(h)x%k 2 ), (z = l,2, j = 2). (2.1) 
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where p\, p 2 and k±, k 2 denote the four momenta of the incoming photons and outgoing 
neutralinos, respectively. In calculation of the amplitude, one should note that there should 
be a relative sign (— l)* 5 ^ between the amplitudes of one diagram and its counterpart obtained 
by exchanging the final neutralinos as a result of Fermi statistics. The corresponding matrix 
element can be written as 

M = M s + M b + M q = e"( Pl )e u (p 2 )u(k 2 )P L {A L 7m7. + tfwrfi + fzlvlv + filulrfi 

+ fnhul^i + fizVi^Piv + fiAP2^Pivi>i + fi^hv + fi^h^i + fi 7 k lfl k lv 
+ fwhpktofa + fi 9 k^p lv + fioh^fa + e^ af3 [f 2 L 1 Pi a p 2 13 + f&Pi a P2%] } v(h) 
+ (p L ^p R ,tf^f*). (2.2) 

where Pl,r = |(1 =F 7s), fk' R {k = 1, ...,22) are form factors. M s ,Mb and M q are the 
amplitudes of the s-channel, box and quartic interaction diagrams, respectively. Since our 
calculation shows that the contribution to the cross section is predominantly from s-channel 
diagrams, we listed the explicit expressions of the form factors for s-channel diagrams in 
Appendix and omitted the form factor expressions of box and quartic interaction diagrams, 
because they have long-winded expressions and negligible contribution to cross section in the 
resonant effect energy region. Then the cross section for this subprocess at one-loop order 
in unpolarized photon collisions can be obtained by 

77 - = ~^— 2 {\P \l dt £|M| 2 , (i = l,j = 1, 2). (2.3) 
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In the above equation, = 1/2 (mi + m% — s) ± \ {m% + m% — s) 2 — 4m| m% 

1 ' ' |_ v Xi Xj 7 y v Xi Xj J Xi Xj 



. The 



factor (!) y is due to the two identical particles in the final states. The bar over the sum 
means average over initial spins. 

The neutralino pair production via photon-photon fusion is only a subprocess of the 
parent e + e~ linear collider. The total cross section of the neutralino pair production via 
photon fusion in e + e~ collider can be obtained by folding the cross section of the subprocess 
(7(77 — > XiXj) with the photon luminosity. 

a(s) = I*™ dz^a( 77 - at s = z 2 s), (2.5) 

where \fs and \fl are the e + e~ and 77 c.m.s. energies respectively and dC 17 /dz is the 
distribution function of photon luminosity, which is 

dC^ fX-max dx , . „ . o 

-jf = 2s / — / 7/e (s)/ 7/e (* 2 /*), (2-6) 



where / 7 / e is the photon structure function of the electron beam [14], [L5|. We take the 



structure function of the photon produced by Compton backscattering as [14, 16 



fComp, \ _ j D )g (l X+ ^ x + £2(^)2 ) , for x < 0.83, 

i7/e K 1 \ 0, for x > 0.83, 



(2.7) 



where 



Taking ujq the maximal energy of backscattering photons, or, £ = 2(1 + y/2), we have 



D(0 = 1.8397. 



3. Numerical results and discussions 



In this section, we present some numerical results of the total cross section in the 
mSUGRA scenario from the complete one-loop diagrams for the process e + e~ — > 77 — > 
XiXj-> (* — 1) j — 1)2). The input parameters are chosen as: m t = 175 GeV, mz = 
91.187 GeV, m b = 4.5 GeV, sin 2 6 W = 0.2315, and a = 1/128. We assume that x°i is the 
LSP and escapes detection. In numerical calculation to get the low energy scenario from 
the mSUGRA, the renormalization group equations(RGE's) |T7|] are run from the weak scale 
Mz up to the GUT scale, taking all thresholds into account. We use two loop RGE's only 
for the gauge couplings and the one-loop RGE's for the other supersymmetric parameters. 
The GUT scale boundary conditions are imposed and the RGE's are run back to Mz, again 
taking threshold into account. 

As we know that the s-channel resonance effects of the intermediate Higgs bosons (see 
Fig.l(b,c)) could play an important role in some c.m.s. energy regions of incoming photons. 
With the mSUGRA scenario input parameters used in our numerical calculation, the mass 
of the lightest Higgs boson h° is under the thresholds of W + W~ and Z°Z° decays. Therefore 
the relevant decay width of h° is mainly contributed by the decay channel of h° — > bb. The 
decay channels of H° may involve H° — ■> gg(where q may be top and bottom quarks) , H° — > 
h°h° (A A ), H° -> xtXi (XiXj(i = 1, j = 1, 2)) and H° -> W+W~ (Z°Z°), if the mass of 
H° is larger than the thresholds of all those decay channels. The main decay channels for the 
pseudo-scalar Higgs boson are A — > Z°h°, A — > qq and A — > XiXi — lj 3 — 1> 2)) 

under similar conditions. In our calculation, all the decay widths of the intermediate Higgs 
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bosons are considered at the tree level and their expressions can be found in appendix B of 



reference ||13|| . 

The cross sections for X1X2 an d X2X2 v * a photon collisions at linear colliders versus 
the mass of x% is shown in Fig. 2. The input parameters are chosen as tuq = 100 GeV, 
A = 300 GeV, /1 > and tan/3 = 4. We calculate the cross sections at the collision 
energies of electron-positron yfs being 500 GeV and 1 TeV, respectively. In framework of 
the mSUGRA, when the mass of \X varies from 81 GeV to 251 GeV as shown in Fig.2, 
chargino mass has the values above the lower limit give in Ref. j?J, the masses of Higgs boson 
H° and A increase from 231 GeV to 510 GeV and from 226 GeV to 507 GeV respectively 
and the mass of the Higgs boson h° is always under the threshold of X1X2 decay. Therefore 
the contribution to the cross section from the s-channel with exchanging h° is very small due 
to the s-channel suppression. The masses of H° and A may reach the thresholds of X1X2 
and X2X2 decay modes, then the cross section of the subprocess will be strongly enhanced 
by the s-channel resonant effects when the total energy of the subprocess approaches the 
mass of H° or A . 

We can see from Fig.2 that all the cross sections for the neutralino pairs x^X^j (i = 
1,2, j = 1) decreases with the increment of the value of m^ 2 . It can reach 18 femto barn 
when y/s = 500 GeV and m^ 2 ~ 81 GeV . For the two curves of \fs = 500 GeV, the cross 
sections with \fs = 500 GeV go down rapidly when the values of the mass of X2 are i n the 
vicinity of 193 GeV. The reason is that the resonant effects from H° and A exchanging 
s-channels do not contribute to the total cross section, since the maximum c.m.s energy of 
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incoming photons is \Js max ~ 0.83^ < m H o(m A o). That fact reflects that the resonant 
effect of the s-channels plays an important role in the total cross sections of neutralino pair 
productions at LC. From this figure we can see that the cross sections of neutralino pair 
productions with y/s = 500 GeV are always larger than those with y/s = 1 TeV when m^o 
is less than 180 GeV, and will be smaller than those with y/~s = 1 TeV when m^ 2 is larger 
than 185 GeV. 

In Fig. 3 we present the cross sections of neutralino pair productions versus tan (3 with the 
colliding energy of electron and positron being 500 GeV and 1 TeV respectively, where the 
other four input parameters are chosen as m = 100 GeV, mi/ 2 = 150 GeV, A = 300 GeV 
and /i > 0. Our calculation in the mSUGRA model shows that when tan/3 increases from 2 
to 32, the mass of x° ranges from 51 GeV to 57 GeV and the mass of xl ranges from 98 GeV 
to 102 GeV. So the dependence of the masses of neutralinos on tan (3 is very weak, but the 
masses of Higgs boson H° and A depend on tan (3 obviously and decrease from 352 GeV 
to 166 GeV and from 344 GeV to 165 GeV , respectively. The mass of h° is a function of 
tan/3 too, but keeps m h o < m^o + m^o. Since the couplings of Higgs bosons to quarks and 
squarks pair of the third generation are related to the ratio of the vacuum expectation values, 
tan (3 should effect the cross sections substantially due to Yukawa coupling strength. The 
line shapes in this figure are determined mainly by parts of the form factors of fi : ' s R , fi^ 
and f which come from the diagrams involving couplings between heavy quarks (top or 
bottom quark) and Higgs bosons. The form factors from these diagrams are all proportional 
to a factor, which has the form as -^-3 + -A3. Therefore all the curves in Figure 3 have 

' cos/3 sm/3 & 
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similar line shapes. We can see from this figure that the cross section for mixed neutralino 
pair (X1X2) production is larger than that for diagonal pairf^XiD production. When the 
ratio of both vacuum expectation values tan j3 is in the vicinity of tan (3 ~ 5, the cross sections 
for both X1X2 an d X2X2 V s * 11 productions can reach maximum, where we have m^ a = 52 GeV 
and m^ 2 = 98 GeV. 

4. Summary 

In this paper, we studied the pair production processes of the neutralinos via photon- 
photon fusion at LC. The numerical analysis of their production rates is carried out in the 
mSUGRA scenario with some typical parameter sets. The results show that the cross section 
of the neutralino pair productions via photon-photon collisions can reach about 18 femto barn 
for X1X2 P a i r production and 9 femto barn for X2X2 P a i r production at an electro-positron 
LC operating in 77 collision mode with our chosen parameters, which is one order smaller 
quantitatively than at the machine operating in e + e~ collision mode. At future TESLA 
collider, the annual 77 luminosity is designed to be 10 — 30 fb~ l and one order higher after 
upgraded[0, these translate into about 90 ~ 270 events per year (about 10 3 events after 
upgrade) for X2X2 production and 180 ~ 540 events per year (few thausand events after 
upgrad) for X1X2 production. Our numerical results present a fact that in some c.m.s energy 
regions of incoming photons, there exist obvious resonance effects from s-channel diagrams by 
exchanging intermediate H° and A Higgs bosons, which can makes observable enhancement 
on the cross sections of the parent processes at linear colliders. 
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Appendix 

In this appendix we list the form factors for the one-loop s-channel diagrams. The 
amplitude of s-channel diagrams, which was defined in Eq.(2.2), has the form 

Ms = -^- i ^{p 1 )e u (p 2 )u(k 2 )P L {/f/vy,, + is L s 7^ + fi3, s P2vPi» + /2 L i, s Wk}"(^i) 
+ (Pl -> Pr, fk,s fk,s) 

e 2 

+ (Pl -> Pr, fk, s ~^ fk, s )i 

where = / 3 L S , /£ = f* s . 

We divide the form factors into four parts: = fk^' W + fkfs^'* + fk^^ + fkfs^'* ■> 

which correspond to loop diagrams of W ± , charged Higgs bosons and ghost particles in 
Fig.l(b,4, b.5 and c.3), chargino loops in Fig.l(c.2), quark and lepton loops in Fig.l(c.2), 
and scalar quark and lepton loops in Fig.l(b.3, c.l) respectively. Notations used in the 
appendix are defined below, 

B mn = P>mn[-Pi,m y , m y ], B% n = B mn [-p 2 , m y , m y ], B% n = B mn \k x + k 2 , m y , my], 
^fnn = C mn [—p2,—pi,m y ,m y ,m y \, 

12 



B 



, / (B = h°,H°,A°), 

S — VTiB + imBL B 



G 



s — m z o* 

The couplings of Higgs boson(B)-Xj-Xj (x = X°, X + ) are denoted as 



where B = h°, H°, A , G°. All the explicit expressions of relevant vertices can be found in 



the Appendix of [13]. The form factors are expressed explicitly as follows, g represents 

electroweak coupling constant; Qf denotes electrical charge(Q T = — 1, Qt = 2/3, etc.) and 
f ( f) 

the factor 3 in f k ^ arises from quark color. 

fxf = Ph {i[V hH+H+ (B™ + - 4C 2 f) + V hG+G+ (Bl w - m w *C™ - 4C%)] 



+ 
+ 



gm w 



L,W 
13,s 



° \ ° + (6 - 3)I?r + (m w 2 + 31)^ - |(C7}T + C 
2C 2 ^) + (16 - 5e)C£] sin (a-0)} V^o 
-> #°, sin (a - /?) -»• - cos (a - /?)), 

8i(-P/i {VfcH+H+(C^ + C23 ) + VhG+G+ (Cl2 + ^23 ) 



w 

12 



+ 
+ 



igm w 



2C™ + 3(C£ + a 



w 



y 23 j 



sin 



(«-/?)Kx°x° 



(h° -> sin (a - /3) -»■ - cos (a - /?)), 



l,x 



ft 



fL,x 
Jl3,s 



12 



24 



x=l 



-4z 2 [cf + 4(C* + + 



23 



2 = 1 



+ -> F°) + -> A ) + (A -> G°) 
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J21,s 



m- 



fi 



LJ 



+ 



+ 



x=l 

,0 , riO^ 



PhV h %o^o(V h L + - - V h R +-) + (h° -> H°) + (h° A ) 



s 



J 13. 



J21,s 



(h° -> G u 

-2«Q r 2 m r (2£? lT - sC T - 2sC[ 2 - 8C T 24 )(P h V hrT V h L ^ + PffV^V^) 
3(r^t)+3(r^6), 

-8*Q r 2 m r [CI + 4(C[ 2 + C 2 r 3 )] (P h V hTT V h L ^o + P h V Htt V^o) 
3(r ^t) + 3(r^6), 

-8m r C£(P A Vk rr V^o + P G K GTr V^o) + 3(r -> t) + 3(r -> 6), 



Jl3,s 



zQ r 2 £(P 3f - ^(Wm^^xO + PHVHw.Vjtqq) + 3(r -> t) + 3(r - 6), 
8*Q T 2 E(C[ 2 + C f 23 )(P h V hfxfx V h L ^ + P H V Hfxfx V^) + 3(r -> t) + 3(r -> 6), 



x=l 



/£ = /£(^x°x° - Vgojo), (i = 1, 13, 215 = 5°, A , G°). 



In the above expressions we adopted the definitions of one-loop integral functions in [18 



and defined d = 4 — 2e. The numerical calculation of the vector and tensor loop integral 



functions can be traced back to four scalar loop integrals A , B , Co, D as shown in [ [19] 
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Figure Captions 

Fig.l The Feynman diagrams of the subprocess 77 — > XiX'j- ( a -l ~ 6) Box diagrams. 
(6.1 ~ 5) Quartic interaction diagrams, (c.l ~ 3) Triangle diagrams. 

Fig. 2 Total cross sections of the process e + e~ — > 77 — > x°x°, = 1? 2, j = 2) as 
function of at electron-positron LC with ^/s = 500 GreV\ The solid curve and dashed 

curve are for X1X2 production at LC with y/s = 500 GeV, 1000 GeV, respectively. The 
dotted curve and dash-dotted curve are for X2X2 production at LC with y/s = 500 GeV and 
1000 GeV, respectively. 

Fig. 3 Total cross sections of the process e + e _ — > 77 — > x°x°, = 1)2, j = 2) as 
function of tan (3. The dashed-lines are for X2X2 production and the full-line curves are for 
X1X2 production, with e + e~ colliding energy \fs = 500 GeV and 1000 GeV respectively. 
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